Spatial Ordering of Defects and Conductivity of Functionalized Graphene 
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Recently, it was predicted that an RKKY-type interaction between adatoms in graphene can 
drive an ordering transition to a state with broken sublattice symmetry ( arXiv: 1004.3678 ). In this 
state, due to Bragg scattering of electron waves on the sublattice modulation, a gap opens up at the 
Dirac point (DP). Here we investigate the effect of ordering on the transport properties, finding that 
upon transition from disordered state to an ordered state the conductivity is increased (reduced) at 
energies away (near) DP. This behavior can be understood as a result of coherent Bragg scattering 
in the ordered state. The conductivity change can serve as a direct signature of adatom ordering. 



I. INTRODUCTION 

Electronic properties of graphene can change in dra- 
matic ways in the presence of adatoms or chemical 
groups, such as H, F or OH, which bind to carbon atoms 
via sp 2 — ¥ sp 3 hybridization, effectively turning the func- 
tionalized carbon atoms into vacancies. Microscopically, 
such defects represent resonant scatterers for the con- 
duction electrons, with the resonance energy positioned 
near the Dirac point (DP) 1 . Resonant enhancement of 
scattering at low energies has dramatic effect on elec- 
tronic properties. In recent experiments on hydrogenated 
graphene 2 ^ it was found that even minute amounts of 
hydrogen can make the resistance near DP very high, 
turning the temperature dependence from metal-like to 
semiconductor-like. This behavior demonstrates that 
adatoms on graphene offer a new knob to control elec- 
tronic properties of this material. 

Properties of adatoms on metal surfaces, and on 
graphene, are strongly affected by the RKKY-type in- 
teractions mediated by the conduction electrons^—. The 
RKKY interaction between adatoms in graphene was pre- 
dicted to be strongly enhanced due to the resonant char- 
acter of scattering at low energies^. In addition, the sign 
of interaction was found to be sublattice dependent, be- 
ing different for the atoms residing on the same or on 
different sublattices, A or B. Such sublattice dependence 
indicates that the RKKY interaction may drive a tran- 
sition to an ordered state in which the adatoms prefer- 
entially populate the A or B sites^^. Bragg scattering 
of electron waves on the sublattice modulation results in 
band gap opening at DP, stabilizing the ordered state by 
a Peierls-type mechanism. The spectral gap was found to 
be nearly free of localized states as long as the adatoms 
reside preferentially on sites of one type, A or B, even 
when their distribution is otherwise random. 

Here we shall discuss how such ordering manifests itself 
in transport properties, focusing on the metallic state re- 
alized far from DP. In our analysis of resonant scattering 
on the adatoms we use the T-matrix approach, similar 
to that used to describe scattering on vacancies^. The 
system with the defects randomly distributed over sub- 
lattices A and B arc treated within the self-consistent 



T-matrix approximation (SCTA), in which we perform 
a virial expansion in the adatom concentration ua, tib- 
The contribution to the conductivity which depends on 
the sublattice occupancy fraction ua/tib arises at sub- 
leading order in the expansion in powers of tia and 
tib, see Eq. pop below. Our virial expansion is valid in 
the metallic regime, at densities which are detuned far 
enough from neutrality, such that 
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with Uf the Fermi momentum and i the mean free path. 

The contribution to conductivity Act which depends 
on ua/tib exhibits different behavior at energies near DP 
and far from it. Near DP we find that ordering results in 
a decrease in conductivity, Act < 0, which is consistent 
with the picture of an energy gap opening at e > near 
DP upon ordering 9 . The contribution Act changes sign 
at an energy e rs 0.5<5, where 5 is an energy parameter de- 
scribing the resonance structure of a T-matrix [see FigQ] 
and Eq.®]. 

In the metallic regime realized at high detuning from 
DP, |e| 3> \6\, we find a positive contribution which ex- 
hibits a universal behavior 



Act( £ ) 



(n A - n B ) 



7T 2 (n A + n B ) 2 \e\ 



(2) 



where W is a bandwidth parameter [see Eq.(|4])]. The pos- 
itive sign of Act and its dependence on ua/tib mean that 
this contribution makes the net conductivity higher at 
n a 7^ tlb than at tia = tib- An increase in conductivity 
due to ordering, Act > 0, may seem somewhat surprizing 
given that the ordered state is characterized by a band 
gap at DP. This behavior can be understood in terms of 
Bragg scattering of electron waves on the disorder poten- 
tial due to the adatoms. In the ordered state, some of 
the harmonics of the disorder potential scatter electrons 
coherently, providing a gap-opening perturbation. Such 
harmonics cease to contribute to the incoherent scatter- 
ing responsible for momentum relaxation. This leads to 
an increase in the mean free path, and in conductivity. 
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II. T-MATRIX OF A SINGLE SCATTERER 

When the adatoms and groups such as H, F and OH 
bind to a carbon atom in graphene lattice, one of the it 
electrons turns into a a electron, thereby being removed 
from the conduction band. As a result, the site of the 
functionalized carbon atom is seen by other it electrons as 
a vacancy. Simultaneously, a quasibound electron state 
apears, giving rise to a low-energy resonance in the lo- 
cal density of states near the adatom. Ab initio study 1 
predicts the resonance energy values which span a wide 
range: = -0.03, -0.11, -0.70, -0.67 eV for H, CH 3 , 
OH and F, respectively. The resonance at e = mani- 
fests itself as a resonance peak in the electron scattering 
amplitude. 

A simplest model for electron scattering on point-like 
defects is provided by a delta-function potential, de- 
scribed by the Hamiltonian 



H= *o(^£V>x'+h.c.)+^w(x)ViV* 



(3) 



|x-x'| = l 



with the hopping amplitude to ~ 3.1 eV, and potential 
w(x) equals U on the site occupied by an adatom and 
zero elsewhere. 

The regime of weak scattering, described by small 
U <C t , can be easily treated in the Born approxima- 
tion (see e.g. Refs. Illfl4l |). Here we are interested in 
the regime of strong scattering, U ^> to- Physically, 
such strong scattering potential corresponds to effectively 
removing an electronic level from the carbon site occu- 
pied by an adatom, by prohibiting electrons to visit it. 
Scattering on a strong, point-like defect can be conve- 
niently described by a T-matrix, defined via the Greens 
function for multiple scattering. In the operator form, 
T = U + TGU , where U is the defect potential and G is 
the free particle Greens function. 

Direct calculation shows that the behavior of the T- 
matrix for the tight-binding model © is the same as that 
found in the continual massless Dirac particle model with 
two degenerate valleys^. Treating the defect as a delta- 
function potential with equal inter- valley and intra- valley 
scattering amplitudes, we find 



T(e) 



An(-W 2 /s 2 



U 



■n 



(4) 



where W 



3t , no 



3V3a 2 



is the density of carbon 



atoms in one sublattice (thus 2no is the net density of 
carbon atoms in graphene lattice). 

The T-matrix Q has a resonant form, with a com- 
plex pole at an energy determined from the equation 
6 + ie* ln(— W 2 /e 2 ) = 0. Since 8 is small for a strong 
scatterer, U 3> to, the energy is approximately given 
by 



The resonance ([5]) is positioned at a low energy, as ex- 
pected from the equivalence between the tight-binding 
and continual models. The energy vanishes, when the 
adatom potential is infinitely strong, U —> oo, turning 
into a zero mode associated with a vacanc y 12 ' 13 . 

In our resonant scattering model, the parameters U 
and S appear as primary quantities which are used to de- 
rive the resonance energy e*. However, since the energy 
value e* is provided by ab initio calculations, we use its 
relation with the parameters S and U to fix their values 
in order to generate a correct value of e*. 

While these results seem reasonable, one may recognize 
that potential strength exceeding the bandwidth param- 
eter to is unphysical, and seek an alternative description. 
In Refs.[l|[l|| the defects arising as a result of sp 2 — s- sp 3 
hybridization were modelled using a localized level hy- 
bridized with the 7r-band of graphene. One can see by 
a direct calculation that the T-matrix in this case has 
the same behavior at low energies as that found for a 
point-like scatterer, Eq.(|l|). 

The agreement between the two models is by no means 
a coincidence. On very general grounds, we can start 
from the unitarity theorem for the T-matrix of a point- 
like scatterer, 



T*(e)-To(e)= i 7T±±\To(e)f 



1TVZ 



- 



(G) 



and 



Rewriting this relation as T (e) — T Q 1 (e) = lit- „ 

taking into account that To(e) is analytic in the upper 
halfplane Ime > 0, we find 



Us) 



2ttv 2 o 



(7) 



where W is the band width, and the function / must 
be purely real for the T-matrix to satisfy the unitarity 
relation. 

In our calculation, Tq(e) will enter the retarded Greens 
function. Accordingly, in the above equation we choose 
the branch of the logarithm that is analytic in the upper 
half-plane, and is real on the imaginary axis. The T- 
matrix which is analytic in the lower half-plane, obtained 
by complex conjugation, is given by Tq (e). This quantity 
will enter the advanced Greens function. 

The above form of the T-matrix, Eq.([7]), is quite gen- 
eral. It includes the T-matrices obtained in the two mod- 
els proposed for description of adatoms - the local poten- 
tial model, Eq.([3]), and the model of a localized level hy- 
bridized with the 7r-band of graphene 1,16 . The local po- 
tential model, described above, corresponds to f(e) = 26, 
where S is defined in Eq.Q. 



III. THE SELF-CONSISTENT T-MATRIX 
APPROXIMATION 



ln(W/|<5|) -ifsgni" 



(5) 



The T-matrix approach can be used to describe elec- 
tron scattering in the presence of an ensemble of adatoms. 
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Generalizing our discussion of a single scatterer, we 
model the system by the Hamiltonian © with the po- 
tential 



t(x) = J2 U5 ( x - x «) 



(8) 



which equals U on the sites occupied by adatoms, and 
zero elsewhere. We shall assume that the adatoms pop- 
ulate the sublattices A and B unequally, with concen- 
trations jia and ub-, while the distribution within each 
sublattice is completely random (see Ref.[9|). The con- 
centrations riA and ns are taken to be equal in the dis- 
ordered state and are unequal in the ordered state. 

To make connection with previous work, we note that 
the dimensionless occupancies < tia < 1 and < ub < 
1 used in Ref.0 correspond to ua/tlq and ris/rio in the 
present notation. The dimcnsionfull concentrations are 
denoted tia and tib in Ref. @. 

Within the SCTA approach, the collective effect of 
disorder on electron scattering can be described by 
the disorder-averaged Greens functions and self-energies 
which satsify a selfconsistency relationii. In our case, the 
disorder-averaged Greens function for each valley takes 
the following form 



G(e,k) = 



1 



Uocrk 



s = 





S B 



(9) 



where a = (cr x ,a y ) are the Pauli matrices. 

The self-energies are obtained by averaging the T- 
matrices over the position of impurities, separately for 
the A and B sublattices, 



^A^n A (T A ), Ss = n B (T E 



(10) 



By a direct calculation we express the disorder-averaged 
T- matrices in Eq. (fT0|) through the self-energies: 



(Ta) 



(Tb 



/(*) + lfl»(-S)' 



/( £ ) + ¥ln(-^)' 
where the quantities ea, £b are defined as 

£a = £ — ^a, Eb = £ — S B . 



(11) 



(12) 



(13) 



The branch of the logarithm in Eqs. pT|) . (|12p is defined 
in the same way as in the definiton of the T-matrix for a 
single scatterer, Eq.®. 

The analogues of Eqs. (HUJ) , (HU , (H2) for the advanced 
Greens function can be obtained by replacing the self- 
energies ^a(b) and the quantities £a(b) by their complex 
conjugates, T,* A{B) and e* A(B) =e-Y,* A{B) . 

The equations for self-energies and can be solved nu- 
merically, and, in the limit of high density, analyti- 
cally (see next section) . The disorder-averaged electronic 
properties of the system, such as density of states and 
conductivity, can be determined from quantities sa,£b- 



IV. CALCULATION OF CONDUCTIVITY 

We seak to evaluate the conductivity of graphene as a 
function of the adatom concentrations on the two sub- 
lattices, tia and jib- We will use a general expression 
for conductivity, obtained from SCTA, in which we per- 
form virial expansion in powers of ua and ns- The lead- 
ing (Drude) contribution is inversely proportional to the 
total concentration ha + tlb, and is thus insensitive to 
the imbalance of ua and ns- In subleading order, we 
find a contribution to conductivity which depends on the 
sublattice imbalance as Act oc uaUb / {tia + ns) 2 - This 
contribution is negative, which makes the conductivity 
of the ordered state, ua ^ tlb, greater than that of the 
disordered state, ua = tlb- 

We start with a general expression for conductivity^, 



a(e) = -e 2 v 2 U RA (e), 

7T 



where e is the Fermi energy, and 

U RA (e) = Tr J ^L G R (e,k)G A (e,k). 



(14) 



(15) 



Here the integration is performed over the whole Brillouin 
zone, and trace is taken over sublattice indices. 

Using the disorder-averaged Greens function (j9|), and 
evaluating the integral over k, we rewrite 



Il RA (e) 



— {b {£ A £b) 



2ttv-> £ A £b 



where F = In 



F* = In - 



V £ A e Bj 



F*(ea£b)), 
(16) 

with the log- 



arithm branch chosen to be analytic in the upper (lower) 
half-plane. In particular, for real £ a = £b = £, we have 



F = In 



- i7rsgn(e); we will 



-Mrsgn(£), F* = In ^ 

use this relation below. 

Eqs. (fT4")) , (|16l) give a general answer for conductivity in 
terms of self-energies. Performing virial expansion of Y>a 
and Ss in ua and ns, we write 



= n A (T A ) - n 2 A (T A ) 



d(T A ) 
de A 



-n A n B (TE 



d(T A ) 
de B 



n B (T B ) - n 2 B (T B ) 



d(T B ) 

dSB 



-jiatib (Ta) 



9(Tb) 



(17) 



(18) 



where the T- matrices (Ta) and (Tb) and their derivatives 
on the right hand side are evaluated at tia = tib = 0, as 
appropriate for Taylor expansion of Eqs. (TIT]) . (IT!?]) . 

At leading order, keeping only the terms of first order 
in ua and ns, we find 



nr ( £ ) = 4 



i 



Vq ua 



i\e\ 

-n B T *(£)-To(£)' 



(19) 
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This expression can be simplified using the unitarity re- 
lation, Eq.©, to obtain 



<T (e) = - 



n\T (E)\ 2 (n A + n B )' 



(20) 



This expression, which is inversely proportional to the 
total adatom concentration ua + n B , is nothing but the 
Drude conductivity. Being a function of the total con- 
centration, this contribution is not sensitive to sublattice 
imbalance. 

For energies e sufficiently far from DP, |e| ^ \e*\, we 
can approximate T(e) « iiv 2 , / (ehi(iW / e)) . In this limit, 
Eq.fl3U|) becomes 



e 2 s 2 ln 2 (W 2 /e 2 ) 
2tt 3 v 2 n A + n B 



(21) 



which is the result found in Ref.jlll]. 

The expansion in powers of ua and n B , employed in 
this and next sections, is controlled by the parameter 



n A{B) \dT/de\<^l. 



(22) 



Evaluating the derivative, dT/de = ttvq ln(iW / e) / (S 
e \n(iW / e)) 2 , we can write this condition as 



^|T ( £ )| 2 m^«i 



(23) 



SawL which cannot be obtained from SCTA. However, 
as we argue below, to leading order this contribution is 
independent of ua and n B and thus is inessential. 

We proceed with Taylor expanding Eqs. ([TU| . (ITT1) . ([T2^ 
to second order in ua and ns as outlined above. In order 
to extract the contribution that depends on the sublattice 
imbalance n A /n B , we rearrange the quadratic terms in 
Eq.flUP as follows: 



<{T A ) 



9{Ta) 



x(Ta) 



de, 



f nAn B (T B 
n A n B ( {Ti 



d(T A 



n A (n A + n B ) 
d(T A 



det 



(Ta 



de A 



and similarly in Eq. p8|) . The terms proportional to 
tia ~\~ n B can be dropped since they cannot produce a 
contribution that depends on the imbalance. Evaluating 
the derivatives in sa and s B , we note that the nonzero 
answer arises solely due to the prefactors Ea and of 
the log terms in Eqs. (ITl"1) .(fl"2" j) . since the contributions 
arising from differentiating ln(— W 2 /sa£b) cancel each 
other. Keeping only the terms proportional to tiati B , we 
write 

Sa(b) = n A{B) T (e) + ^ln A n B F(e), (24) 
^A(B) = n A{B) n{e) + Tt^-n A n B F*(e), (25) 



This condition can also be written in terms of the mean 
free path, see Eq.([T]). Recalling that Drude conductiv- 
ity in two dimensions is of order kp£e 2 /h, and com- 
paring to the result (|2T)]) . we obtain Eq.([T|). Thus, no- 
tably, the parameter controlling our virial expansion is 
k F £/\n(W/\e\) > I rather than k F l > 1. 



V. THE CONTRIBUTION Act 

The subleading terms in the expansion of conductiv- 
ity in powers of adatom concentration are zero-order in 
ua and ns- According to Eq. ([!]), such terms are also 
zero-order in k F £ and are of order e 2 /h. As we shall see, 
there are two distinct contributions to conductivitiy aris- 
ing at this order. One is a contribution proportional to 
tiait-b I '(n A + ns) 2 , which we shall denote as Act. This 
contribution depends on the imbalance n A /n,B and thus 
describes the effect of sublattice ordering on conductiv- 
ity. In addition, there is a weak localization contribution 



These expressions are then plugged into Eq. (fl6|) , and ex- 
panded to keep the two lowest orders in n A , ub- This 
gives 



T RA 



HA i 



(26) 



where H RA (e) is given by Eq. (fT9|) . The contribution H RA 
can be represented as a sum of two terms, 



nf A ( £ ) = nf A ( e ) + (5nf A (£ 



(27) 



where II ^ depends on n A and while JII^" 4 does not. 
The latter term, which is given by 

ttv 2 v 2 T*( £ )-T ( £ )' V ' 

will be ignored in what follows. 

The contribution fl^ , after some algebra, can be writ- 
ten as 



RA 



2i\e\ 



v 2 (T*(e)-T (e)) 2 (n A 



UATlB 



- n B ) 



T 2 (e)-T* 2 (e) 



«sgne 3 
— 2-( T o ( £ ) 



The three terms in this expression originate from an expansion of the denominator in Eq. (|I6p in n A , ns 
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as follows: the first term comes from the products in 
S^Es — T,* A T,* B , where self-energies are taken in the low- 
est order in tia, tib', the second and third terms come 
from the term — e(IU + Eb — — E B ), where the self- 
energies are taken in the second order in tia, tib- 



The result (|29|) can be simplified using the unitarity 
relation for the T- matrix, Eq.(j6|). Subtracting the value 
at ha — ub, and using Eq.([T4|). gives a contribution 
describing the change in conductivity due to ordering, 



Act 



e 2 vl (n A - n B ) 2 
27r|T (e)| 2 (n A +n B ) 2 



e\To{e) 



m 



Energy dependence of Act is shown in Fig[TJ The contri- 
bution (|30|) is positive sufficiently far from DP, |e| ^> |e* [, 
and negative near DP at e > 0. 

The expression (|30l) simplifies at large detuning from 
DP, since at large e the imaginary part of Tq(s) is small 
compared to its real part. Thus, in Eq. (|30|) we can ap- 
proximate T (e) w T *(e) » ReT (e). This gives 



A<7 



2 2 



(n A -~n B ) 2 



27r(ReT ( £ )) 2 (n A + n B ) 2 



(31) 



-4 



ReT (e) 



— (ReT (e)YHW z /e z ) 

?™0 



Replacing ReT (e) by nvQ/(e\n(W/\e\)), we obtain the 
desired contribution to the conductivity, Eq.Q. 

Taking the ratio of Act and the Drude conductivity, 
Eas.([2"]).([2"T]). we have 




-1 o 1 

Energy e/8 



Act HW/\e\) 
ctq kpi 



(32) 



where we used the relation oq ~ fcp-f e 2 //i. Comparing to 
Eq.([T]), we see that Act/cto <C 1 in the regime of validity 
of our virial expansion. 

It is interesting to compare the contribution Act, 
Eq.([2]), to the weak localization (WL) contribution, 
which arises in the same order of an expansion in 1/ (hpl). 
The WL contribution is given byii 



Sctwl = --^ln— , 



(33) 



where L is the sample size. We note that at leading order 
in 1 / (kp€) the mean free path I depends only on the total 
adatom concentration ua + tib , but not on the sublattice 
imbalance. This can be understood by combining the 
result for Drude conductivity, Eq. (|2"T1) . with the relation 
kp£ = crh/e 2 . We obtain, to leading order in /cf£, 



FIG. 1: Energy dependence of the contribution to conductiv- 
ity which is sensitive to sublattice ordering. Shown is the con- 
tribution obtained by virial expansion in ua, tib, Eq. (|3U|l . and 
its high-energy asymptotic form, Eq.©, for W/5 = 20. Real 
and imaginary parts of a single-scatterer T-matrix, Eq.@, 
are shown by dashed lines. The virial expansion is valid at 
densities detuned far enough from DP, Eq.JTJ). Conductivity 
increases due to sublattice ordering at energies away from DP 
and is reduced at energies e > near DP, in the region where 
an energy gap opens up in the ordered state. 

imbalance ua/tib- However, the corresponding contribu- 
tion to the 8owl, Eq. ([33|) . arises only from I under the 
log, and thus is of order l/(kp£) in units e 2 /h. This is 
parametrically smaller than the contribution Act, Eq.@. 
From this we conclude that the effect of sublattice im- 
balance arising from the WL contribution is inessential 
compared to the Act contribution. 



VI. SUMMARY AND DISCUSSION 



to 



Tr 2 v (n A + n B ) 



ln 2 (W 2 /s 2 ), 



(34) 



a result which is independent of the sublattice imbalance. 

The subleading term in the expansion in powers of kp£, 
which is of order 6£ ~ £j{kp€) = may depend on the 



In summary, sublattice ordering of defects directly 
manifests itself in transport properties. Conductivity in- 
creases in the metallic regime far from DP, because in 
the ordered state part of the scattering ceases to con- 
tribute to momentum relaxation, instead acquiring a co- 
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herent, Bragg-like character. Simultaneously, in the or- 
dered state, conductivity is reduced at energies near DP, 
where the band gap opens up£. 

These effects were investigated using a virial expansion 
of conductivity in the adatom concentration, a method 
that is valid in the metallic state at sufficiently large 
mean free path, Eq.([T]). One can expect transport in 
the state with a well developed gap to be dominated by 
the low-energy modes propagating along domain bound- 
aries, separating regions with opposite signs of the order 
parameter (sublattice imbalance), as well as by thermal 
activation across the gap. Theoretical description of this 
regime lies outside the scope of this paper. 

Finally, we note that the results of this work apply to 
other scenarios of partial ordering of defects, in which 
populations on several sublattices become spontaneously 
imbalanced, the distribution on each sublattice remaining 



random. One example is the recently discussed Kekule- 
type ordering^ in which three equivalent sublattices be- 
come unequally populated. We expect that the picture of 
Bragg scattering on the density harmonics arising due to 
ordering, accompanied by a reduction in disorder scatter- 
ing and an increase in conductivity, applies in this case. 
An increase in conductivity can serve as an experimental 
signature of a transition to an ordered state. 
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